Novel shell-model analysis of the 136 Xe double beta decay nuclear matrix elements 
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Neutrinoless double beta decay, if observed, could distinguish whether neutrino is a Dirac or a 
Majorana particle, and it could be used to determine the absolute scale of the neutrino masses. 
136 Xe is one of the most promising candidates for observing this rare event. However, until recently 
there were no positive result for the allowed and less rare two-neutrino double beta decay mode. The 
small nuclear matrix element associated with the small half-life represents a challenge for nuclear 
structure models used for its calculation. We report a new shell-model analysis of the two-neutrino 
double beta decay of 136 Xe, which takes into account all relevant nuclear orbitals necessary to fully 
describe the associated Gamow- Teller strength. We further use the new model to analyze the main 
contributions to the neutrinoless double beta decay half-life, and show that they are also diminished. 

PACS numbers: 23.40.Bw, 21.60.Cs, 23.40.Hc 



Neutrinoless double beta (Of/3/3) decay can only occur 
by violating the conservation of the total lepton number, 
and if observed it would unravel physics beyond Standard 
Model (SM) of particle physics and it would represent a 
major milestone in the study of the fundamental proper- 
ties of neutrinos [l| . Recent results from neutrino oscilla- 
tion experiments have demonstrated that neutrinos have 
mass and they mix In addition, they show that 

neutrinoless double beta decay process could be used to 
determine the absolute scale of the neutrino masses, and 
can distinguish whether neutrino is a Dirac or a Majo- 
rana particle [f| . A key ingredient for extracting the ab- 
solute neutrino masses from Of/3/3 decay experiments is a 
precise knowledge of the nuclear matrix elements (NME) 
for this process. There is a large experimental effort in 
US and worldwide to investigate the double beta decay 
of some even-even nuclei [H . Experimental data for two- 
neutrino double-beta decay (2f/3/3) to the ground state 
(g.s.) and excited states already exist for a group of nu- 
clei [6j. There is no confirmed experimental data so far 
for neutrinoless double-beta decay. The prediction, anal- 
ysis and interpretation of experimental results, present 
and expected, are very much dependent on the precise 
nuclear structure calculations of corresponding transition 
probabilities. 

Although many experimental efforts in US and world- 
wide, such as MAJORANA and GERDA Q, are pin- 
pointing to the /3/3 decay of 76 Ge there are very en- 
couraging results related to the /3/3 decay of 136 Xe. For 
a long time there were available only upper limits for 
the 2f/3/3 half-life. Recently, the EXO-200 collabora- 
tion reported @, @] a precise measurement of this half 
life of 2.11 ± 0M{stat) ± 0.21(sys) x 10 21 yr, corre- 
sponding to a NME of 0.019 MeV" 1 . This large half- 



life would imply a smaller background for the associ- 
ated Of/3/3 measurement and EXO, a larger version of 
EXO-200 designed for reaching this goal, is under con- 
sideration [l|. The upper limit for the Of/3/3 half-life re- 
ported by EXO-200 is 1.6 x 10 25 yr @]. In addition, the 
KamLAND-Zen collaboration reported a 2f /3/3 half-life of 
2.38 ± 0.02(stat) ± O.U(sys) x 10 21 yr and a upper-limit 
for the Of /3j3 half-life of 5.7 x 10 24 yr [§]. 

Since most of the /3/3 decay emitters are open shell nu- 
clei, many calculations of the NME have been performed 
within the pnQRPA approach and its extensions U 13 ■ 
However, the pnQRPA calculations of the more common 
two-neutrino double beta decay observable, which were 
measured for about 10 cases [6j, are very sensitive to 
the variation of the g pp parameter (the strength of the 
particle-particle interactions in the 1 + channel) [13l. Il4j. 
and this drawback persists in spite of various improve- 
ments brought by its extensions, including higher-order 
QRPA approaches [12j . Although the QRPA methods 
do not seem to be suited to predict the 2f /3/3 decay half- 
lives, they can use the measured 2f /3/3 decay half-lives to 
calibrate the g pp parameters that are further used to cal- 

Another method that 
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culate the Of /3/3 decay NME 
was recently used to calculate NMEs for most Of /3/3 decay 
cases of interest is the Interactive Boson Model (IBA-2) 
[loT ]. However, a reliable IBA-2 approach for 2f/3/3 decay 
is not yet available. 

Recent progress in computer power, numerical algo- 
rithms, and improved nucleon-nucleon effective interac- 
tions, made possible large-scale configuration-interaction 
(CI) calculations (also known as shell-model calculations) 



of the 2f/3/3 [T3-[3 and Of/3/3 decay NME [H E[|. The 
main advantage of the large-scale shell-model calcula- 
tions is that they take into account all of the many- 
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body correlations for the orbitals near the Fermi sur- 
face. Also they are also less dependent on the effective 
interaction used, as long as these are based on realistic 
nucleon-nucleon interactions with minimal adjustments 
to the single-particle energies and some two-body matrix 
elements so they reproduce general spectroscopy of the 
nuclei involved in the decay [2l| . Their main drawback is 
the limitation imposed by the exploding CI dimensions 
even for limited increase in the size of the valence space 
used. The most important success of the large-scale shell- 
model calculations was the correct prediction of the 2v(3(3 
decay half-life for 48 Ca nam. In addition, the CI calcu- 
lations do not have to adjust any additional parameters, 
i.e. given the effective interaction and the Gamow- Teller 
(GT) quenching factor extracted from the overall spec- 
troscopy in the respective mass-region, they are able to 
reliably predict the 2v(3(3 decay half-life of 48 Ca. 

CI methods provide realistic many-body wave func- 
tions (w.f.) for many nuclei from 16 to 100 Sn and be- 
yond. These wave functions can describe observables 
related to specific experiments, e.g. for nuclear astro- 
physics, and the electro- weak interactions with the nu- 
cleus. The minimal valence space required for 136 Xe 
involves the 0g7/2ld5/2lrf3/22si/20/iii/2 orbitals for pro- 
tons and neutrons (the jj55 model space). There are no 
spurious center-of-mass (CoM) states in the jj'55 model 
space since the CoM operator R does not connect any 
of the orbitals. The key is to obtain effective interac- 
tions (EI) that can provide energies and wave functions 
in jj'55 model space that are at a similar level of accuracy 
as those obtained for the sd-shell [23| and for the p/-shell 

24 1. The CI (3(3 decay NME were reported over the years 
considering continuous improvements of the 
EI. These calculations indicate a significant sensitivity of 
the results to the improving EI. For example, the quench- 
ing factors used to describe 2v(3(3 NME varies from 0.74 

13 to 0.45 Hi], and the 0^/3/3 NME varies by a fac- 
tor of about 3 between Ref. [17| and the more recent 
Ref. [20(| . One of the drawbacks of model spaces such 
as jj'55 is that in order to maintain center-of-mass pu- 
rity they do not include all spin-orbit partners of orbitals 
such as 0g 7 /2 and 0/i n / 2 . The known effect is that the 
Ikeda sum-rule is not satisfied indicating that some the 
Gamow- Teller strength, which is so important for both 
types of NME, is missing from this model space. For 
example, in jj'55 typical Ikeda sum-rule for 136 Xe is 52, 
while the expected result is 84 (see also Table Q] below). 

In this letter we investigate the effect of extending the 
model space to jj77 by including the effects of the miss- 
ing Ogg/2 and O/ig/2 orbitals. The two-body matrix ele- 
ments with good J and T were obtained from the code 
CENS [26]. The procedure discussed below was used to 
obtain a Hamiltonian for the jj77 model space that we 
will refer to as jj 77 a. In the first step, the short-range 
part of the N 3 LO pot ential [27J was integrated out using 
the Viowfe method [28j . The relative two-body matrix ele- 



ments were evaluated in a harmonic-oscillator basis with 
?iu;=7.874 (a value appropriate for 132 Sn). In the second 
step the interaction was renormalized into the jj77 model 
space assuming a 100 Sn closed core. The 0g g / 2 orbital 
was treated as a hole state, while the other are treated 
as particle states. For the energy denominators we take 
all orbits in the jj77 space to be degenerate with the 
other orbitals spaced in units of Tiu above and below. The 
core-polarization calculation used the Q-box method and 
includes all non-folded diagrams through second-order in 
the interaction and sums up the folded diagrams to infi- 
nite order [2{| . Particle-hole excitations up through ATiu 
were included. Matrix elements obtained in the proton- 
neutron basis were transformed to a good-T basis by us- 
ing the neutron-neutron matrix elements for the T = 1 
components. 

The single-particle matrix elements were obtained 
starting with the jj55 model space for a 132 Sn closed 
core. The five single-particle energies for 0^7/2, 1^5/2, 
lc/3/2, 2si/2 and 0hn/ 2 were adjusted to reproduce the 
experimental values for neutron holes related to the spec- 
turm of 131 Sn as given in [3(|. The results obtained for 
the single-particle energies of protons related to the spec- 
trum of 133 Sb are in reasonable agreement with experi- 
ment [3(3] except that the W5/2 energy is too high by 1.2 
MeV and the l/in/2 energy is too high by 2.4 MeV. Re- 
duction of the diagonal two-body matrix elements by 0.3 
MeV for these two orbitals improves the agreement with 
experiment with minimal overall change to the Hamilto- 
nian. The adjustment of the single-particle energies to 
experiment implicitly includes most of the effects due to 
three-body interactions. 

The two-hole spectrum for 130 Sn and the two-particle 
spectrum for 134 Te are in best overall agreement with 
experiment if the T = 1 matrix elements are multiplied 
by 0.9. The results (experiment vs theory) are (1.28, 
1.34) MeV for 130 Sn and (1.22, 1.35) MeV for 134 Te. For 
application to the larger jj77 model space the single- 
neutron hole energy for Ogg / 2 was placed six MeV below 
the 0<?7/2 energy in 131 Sn, and the single-proton particle 
energy for 0h 9 / 2 was placed six MeV above the 0hu/ 2 
energy in 133 Sb. 

The 2i/(3(3 half-life for the transition from the + gjs. 
of 136 Xe to the 0+ g.s. of 136 Ba can be calculated [31 1 
using 
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where G 2v is a phase space factor that for the the 2v(3(3 
of 136 Xe is 1.279 x 10~ 18 yr^MeV 2 (3l|, M^(0+) is the 
2v(3(3 matrix element given by the double Gamow- Teller 
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TABLE I: Matrix elements in MeV -1 for 2v decay calculated 
using the standard quenching factor 0.74 for the Gamow- 
Teller operator using different number of excitations from 
jj55 to the larger model space. Last column indicate the 
calculated Ikeda sum-rule for 136 Xe. 
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0.062 
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0.091 


84 
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0.037 


84 


1 
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0.020 


84 



Here Ek is the excitation energy of the it state of 136 Cs 
and E = iQ^(0+) + AM = 1.31 MeV, where we used 
the recently reported [13] Q- value Qpp(0 + ) = 2.458 MeV 
corresponding to the /3/3 decays to the g.s. of 136 Ba; AM 
is the 136 Cs - 136 Xe mass difference. 

In Ref. [r| we fully diagonalized 250 1 + states in the 
intermediate nucleus to calculate the 2^/3/3 decay NME 
for 48 Ca. This procedure can be used for somewhat heav- 
ier nuclei using the J-scheme shell-model code NuShellX 
[35 | . but for cases with large dimension one needs an al- 
ternative method. Here we used a novel improvement 
[34j ] of the known strength- function approach [16| , which 
is very efficient for large model cases, such as jj'55 and 
jj77. For example, to calculate the NME for the decays 
of 128 Te in jj55 and 136 Xe in jj77 {n = 1 for 0+ and 
n = 2 for 1 + in Table HJ one needs to solve problems 
with m-scheme dimensions of up the order of to ten bil- 
lion. 

The result when restricting the jj77 model space to 
jjhh is given on the first line in Table HI As already men- 
tioned, the Ikeda sum-rule is only 52 rather then 84, 
indicating that not all GT strength is available in the 
jj"55 space. Although the excitation energies of the GT 
strength distribution are reasonably well reproduced, the 
GT operator err has to be multiplied by a quenching fac- 
tor due to correlations beyond the jj77 model space. In 
typical one major-shell calculations, such as the sd or pf, 
this quenching factor was determined to be around 0.74- 
0.77 (see e.g. Ref. [35j, |36j) that is consistent with that 
obtained in second-order perturbation theory (3^, l38j |. 
Here we use 0.74. Ref. [25| suggests that one should 
use a lower quenching factor in the j'j55 model space, 
0.45, to get an NME consistent with the recent exper- 
imental data. Indeed, our matrix elements in the j'j'55 
model space becomes 0.022 MeV -1 when 0.45 is used. 

However, it would be important to check if the missing 
spin-orbit partners are responsible for the larger result; 
the relative phases in Eq. ([2]) could lead to large cancel- 
lations. Here we consider the larger jj77 model space, 
but we could only allowed few n particle being excited 



from the Ogg/2 orbital or to the h g / 2 orbital, relative to 
Table U also presents the NME for different combi- 

TABLE II: Matrix elements for 0v decay using two SRC mod- 
els CD-Bonn (SRC1) and Argonne (SRC2). The upper 
values of the neutrino physics parameters r]j P in units of 10 -7 
are calculated using the G°" from Refs. [U and 0. 
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SRC1 


2.21 


143.0 


1106. 


206.8 




SRC2 


2.06 


98.79 


849.0 


197.2 


n = 1 


SRC1 


1.46 


128.0 


1007 


157.8 




\«7\ Oil 


8.19 


0.093 


0.012 


0.075 




\v7\ [43] 


9.02 


0.103 


0.013 


0.083 



nations of the allowed n for the initial and final + states 
and the intermediate 1 + states. One can see that when 
n is 1 for the + states and 2 for the 1 + states the NME 
decreases almost to the experimental value without the 
need of artificially reducing the quenching factor. In ad- 
dition, the Ikeda sum-rule is always satisfied in the larger 
model space. 

One should mention that in the jj77 model space the 
wave functions could have CoM spurious components. 
We checked our initial and final + g.s. w.f. and we 
found negligible (less than 3 keV) spurious contribution 
to expectation values of the CoM Hamiltonian. We did 
not check the amount of CoM spuriously in the interme- 
diate 1 + states, but it's unlikely to be large because the 
strength function method [lH performs a small number 
of Lanczos iterations (about 30) starting with a door- 
way state obtained by applying the GT operator on the 
largely nonspuroius + state. As a further check we com- 
pared the GT strength (BGT) for the transition from the 
g.s. of 136 Xe to the first 1 + state in 136 Cs with recent 
experimental data |39j. Table I of Ref. [3j| provides a 
BGT of 0.149(21) for the first 1+ state at 0.59 MeV, but 
we learned [40j that this will be updated to 0.24(7). Our 
BGT is 0.51 in the jj55 model space, but 0.34 in the 
largest jj77 model space, much closer to the experimen- 
tal value. Although, we cannot verify if the calculations 
are converged we can conclude that including all spin- 
orbit partners is essential for a good description of the 
2v/3/3 for 136 Xe. 

Having tuned our nuclear structure techniques to get- 
ting an accurate description of the two-neutrino double- 
beta decay we calculate the NME necessary for the anal- 
ysis of the neutrinoless double-beta decay half-life 136 Xe 



2 it , |4l|. Considering the most important mechanisms 
that could be responsible for 0^/3/3 decay [13] one can 
write the 0^/3/3 half-life 
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where M® v NME and r/j neutrino physics parameters for 
light neutrino exchange (j — v) , heavy neutrino exchange 
(j = N), gluino exchange (j = A') and squark- neutrino 
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mechanism (j = N) are described in Refs. 
is a phase space factor tabulated in several publications. 
One widely used value [lH is 43.7 x 1CP 15 yr^ 1 . A recent 
publication [43| proposes 36.05 x 1CP 15 yr _1 , which is 
about 20% lower. The results for the NME calculated in 
the closure approximation are presented in Table HI1 using 
the n — and n = 1 + w.f. (see Table I). Two recent 
short-range correlations (SRC) parametrizations are used 
12 . \2l\ . No quenching of the bare transition operator 



21|, |44|. The My for the jj55 model space 



was used 

(n = 0) is consistent with other recent shell-model results 
[20j . The NME for the other three mechanisms calculated 
within a shell-model approach are reported here for the 
first time. The NME in the largest space (n = 1) are 
10-30% lower. These results suggest that the inclusion of 
the spin-orbit partners, which proved to be significant for 
a good description of the 2v(3j3 NME, could be also very 
important for a reliable description of the 0z//3/3 NME. In 
addition, they indicate that the net effect is a decrease of 
the NME rather than an increase (an assumption often 
used to understand the lower shell-model value relative 
to the results of other methods, such as QRPA, IBA-2, 
Projected Hatree-Fock Bogoliubov [45j, and Generator 
Coordinate Method [3]). Table HT1 also presents upper 
limits for the neutrino physics parameters | iy p | under the 
assumption of single mechanism dominance. They were 
obtained from Eq. © using the lower limit for the half- 
life 1.6 x 10 25 yr from Ref. Q and the two phase space 
factors of Refs. 31 1 and [43|. Using the upper limits for 
\Vvl\ = m i3p/ m e one can extract an upper limit for the 
effective neutrino mass mpp of 0.42-0.46 eV. 



In conclusion, we reported a new shell-model analysis 
of the two-neutrino double beta decay of 136 Xe that takes 
into account all relevant nuclear orbitals necessary for a 
good description of the Gamow- Teller strength. We show 
that this extension of the valence space can account for 
the small NME without recourse to an artificially small 
quenching factor. We also show that it could lead to 
smaller NME for the most interesting neutrinoless double 
beta decay mode. 
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